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A strong-to- weak-coupling duality is established for the nonequilibrium interacting resonant-level 
model, describing tunneling through a single spinless level, capacitively coupled to two leads by a 
contact interaction. For large capacitive coupling, U, and in the presence of a finite voltage bias, 
the model maps onto an equivalent model with a small capacitive coupling pU' ~ ^/pU. Here p is 
the conduction-electron density of states. This duality is generic to all high-energy cutoff schemes, 
however its details may vary from one regularization scheme to another. In particular, it has the 
status of an exact mapping within Abelian bosonization, applicable to all repulsive interactions 
whether weak, intermediate or strong. On a lattice the mapping is restricted to large U, and is 
controlled by the small parameter 1/pU. Explicit expressions are given for the low-energy scale, 
differential conductance, and occupancy of the level in the limit where U is large. 
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I. INTRODUCTION 

The remarkable progress in nanofabrication techniques 
has focused considerable experimental and theoretical in- 
terest on quantum impurity systems out of thermal equi- 
librium. Typically the impurity is realized by a small 
nanostructure, e.g., a semiconductor quantum dot or 
an individual molecule, which is attached to two leads. 
Application of a voltage bias between the leads causes 
steady-state current to flow through the device. Being 
strongly correlated, such systems cannot be studied per- 
turbatively over the full range of parameters, e.g., tem- 
perature or voltage bias, which calls for the develop- 
ment of nonperturbative techniques capable of treating 
the nonequilibrium steady state. 

In contrast to the plethora of methods available for 
tackling quantum impurities in thermal equilibrium, only 
limited progress has been made out of equilibrium. The 
difficulty lies in the need to work directly with open sys- 
tems, subject to appropriate boundary conditions im- 
posed by the bias. The latter boundary conditions can 
be incorporated either at the level of the many-particle 
wave function through solution of a suitable Lippmann- 
Schwinger equation (the time-independent formulation), 
or by adopting a two-operator description as in the 
Keldysh approach^ (the time-dependent formulation). 
Unfortunately, most of the powerful approaches that 
are available in equilibrium are presently inadequate for 
treating open systems. In particular, leading numeri- 
cal methods such as Wilson's numerical renormalization 
group^ or the density-matrix renormalization group^ are 
presently confined to closed finite-size systems. 

Recently, a step forward was taken by Mehta and An- 
drei;^ who generalized the Bethe ansatz approach to al- 
low the description of nonequilibrium systems in a steady 
state. In the traditional thermodynamic Bethe ansatz 
(TBA), one typically works with closed systems, namely. 



finite-size systems subject to periodic boundary condi- 
tions. The limit L ^ oo, L being the size of the system, 
is taken only at the end of the construction. This ap- 
proach, while yielding all the thermodynamic quantities, 
does not allow for the description of the nonequilibrium 
steady state that develops in a biased two-lead system, 
where particles cross the impurity and dissipate their en- 
ergy asymptotically in the leads. A description of such 
processes requires different asymptotic behaviors in the 
past and in the future, represented by x — > ±00 in the 
time-independent formulation. Thus, in the scattering 
Bethe ansatz (SBA) approach"* one works directly with 
open systems, where the limit L ^ 00 is implemented 
from the outset. The many-particle scattering states de- 
fined on the open system satisfy the Lippmann-Schwinger 
equation with asymptotic boundary conditions of free bi- 
ased leads at x ^ — 00. 

As a first demonstration of the approach, Mehta and 
Andrei applied it to the interacting resonant-level model 
(IRLM))^ describing tunneling through a single spinless 
level with capacitive coupling to the leads. The exact 
SBA solution revealed several surprising features, includ- 
ing a nonmonotonic dependence of the current on the 
capacitive coupling U, and an apparent duality between 
[7 = and U ^ 00. These results have stimulated con- 
siderable interest in the nonequilibrium IRLM,^'^'^'^ pro- 
viding complementary points of view and benchmarks for 
comparison. 

In this paper we provide such a benchmark, by estab- 
lishing a strong-to- weak-coupling duality for the nonequi- 
librium IRLM. In equilibrium, duality means that one 
can map the effective low-energy Hamiltonians for large 
and small U. Out of equilibrium the mapping must also 
preserve the boundary conditions imposed by the bias, 
which sets a far more stringent constraint on the map- 
ping. In the following we show that duality between 
strong and weak coupling is a generic feature of the 
nonequilibrium IRLM, independent of the cutoff scheme 
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used. However, its details may vary from one regular- 
ization scheme to another. This is an important point 
to bear in mind whenever comparing different computa- 
tional schemes. Indeed, within Abelian bosonization we 
find an exact mapping between U and U' = A/{'n'^p^U), 
where p is the conduction-electron density of states per 
lattice site. Thus, C/ = and U oo are equivalent 
within bosonization, in agreement with the SBA. A sim- 
ilar duality between pU 3> 1 and pU' ^ ^/pU ^ 1 is 
found on a lattice, however the mapping is only approx- 
imate, and is accompanied by a significant reduction of 
the associated low-energy scale according to To{pU » 
1) ~ To{U = 0)/{pU)'^. Explicit analytical expressions 
are given in the limit pU 3> 1 for the low-energy scale. 
To, the differential conductance, G, and the occupation 
of the level, Ud- 

The outline of the paper is as follows. In Sec. [Til we 
introduce the noncquilibrium IRLM, and specify our for- 
mulation of the finite bias. Using Abelian bosonization, 
an exact mapping between U and U' — 4/(7r^p^[/) is 
derived in Sec. IIIIl followed by treatment of a general 
lattice cutoff in Sec. IIVI Explicit analytic expressions 
are obtained in turn in Sec. |V]for the low-energy scale, 
differential conductance, and occupancy of the level, in 
the limit where U is large. Finally, we present our con- 
clusions in Sec. IVII 



II. 



Hamiltonian 



NONEQUILIBRIUM INTERACTING 
RESONANT-LEVEL MODEL 



As already indicated above, there are several possi- 
ble ways to formulate quantum impurity systems out of 
equilibrium. Here we adopt a two-operator formalism, 
consisting of the Hamiltonian of the unbiased system, 
Tl, and the noncquilibrium operator, Yq, which specifies 
the bias. This choice of operators arises naturally when 
starting from a coupled system in thermal equilibrium, 
and switching on the bias at some remote time Iq. Thus, 
the system is described at time Iq by the density op- 
erator po — e~''^/Tr{e~'^^}, and evolves henceforth in 
time according to the Hamiltonian of the biased system, 
TC = TC+Yq. Steady state is reached if the limit to ^ —oo 
exists, as can be shown in this case.— 

An equivalent, more common formulation of steady 
state starts from two decoupled leads, each in thermal 
equilibrium with its own chemical potential. Tunneling 
is then switched on at time to, which drives the system to 
a new steady state. Both formulations illustrated above 
give identical resultsi^ However, we shall focus on the 
former one, as it affords a more concise presentation. 
In the SBA approach, applied directly in the steady- 
state limit, the second formulation was adopted^ with 
Pi„iti^i = e-/5(Wo-Yo)/Tr{e-''(^o-'>^o)}. Hereto describes 
the biased system with the leads decoupled. 

With the former formulation in mind, the noncqui- 
librium IRLM is defined in the continuum limit by the 



il;j{x)dx'4'j{x)dx 
+ €ddU+ Va^i {V'](0)d + H.c.| 

] = 1.2 

+ UafdU-^) Yl ^V'] (0)^,(0): 

^ ^ 3 = 1,2 



and noncquilibrium operator 
eV 



Yn = 



dx. 



(1) 



(2) 



Here, the left-moving fields with j = 1,2 describe 

free conduction electrons in lead j; creates an electron 
on the level; a is a short-distance cutoff corresponding 
to a lattice spacing; tj is the tunneling matrix element 
between the level and lead j; U is the capacitive coupling 
(contact interaction) between the level and the leads; and 
V is the applied voltage bias. Note that we allow for 
different tunneling matrix elements to each of the leads, 
but the capacitive coupling is taken to be equal. As for 
the fields tj)^^{x), these obey canonical anticommutation 
relations 



{il}dx),'ip]{y)} = 5ij5{x - y), 



(3) 



subject to the regularization 5(0) = 1/a. Hence, the 
conversion between a Wilson-type tight-binding repre- 
sentation^ and the continuum limit is specified by the 
relation cjg = ■y/aV'j(0), where c] q creates a localized 
electron (zeroth Wilson shell) on lead j. As for the sym- 
bol '.Tpjipj :, it stands for normal ordering with respect to 
the unperturbed Fermi seas of the two leads. 



III. ABELIAN BOSONIZATION 

We first treat the nonequilibrium problem defined in 
Eqs. H]) and ([2]) using Abelian bosonization. In the stan- 
dard fashion,^" two bosonic fields ^j{x) are introduced, 
one boson field for each left-moving fermion field. The 
fermion fields are written as 



i^jix) 



-Hj (x) 



V27 



where the $'s obey 



[$j(a;),$j(?/)] ^ -i%7rsign(x- y). 



(4) 



(5) 



The ultraviolet momentum cutoff — n/a is related 
to the conduction-electron bandwidth D and the density 
of states per lattice site p through D = hvp/a and p = 
1/{2D) — a/{2hvp), respectively. The operator e**^ in 
Eq. is a phase-factor operator, which comes to ensure 
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that the different fermions anticommute. Our exphcit 
choice for (p reads 



and 



Ni- N2 + 2d)d 



(6) 



where Nj is the total number operator for electrons in 
lead j. Alternatively, one can replace the phase- factor 
operator appearing in Eq. Q with a Majorana fermion. 

In terms of the boson fields, the Hamiltonian Ti and 
nonequilibrium operator take the forms 



n 



+ 



flVF 

E 



/oo 
-00 

^{e'*^(")e-*^d + H.c.} 



E 



and 



where 



eV 
47r 



[V$i(x)- V$2(x)] 



fTTpU\ 

arctan 

V 2 ; 



(7) 



(8) 



(9) 



is the phase shift associated with U. This bosonized form 
of the U interaction term stems from the cutoff scheme 
used in bosonization."'^^ It follows from the requirement 
that the exact scattering phase shift be reproduced in the 
limit ti,t2 — > 0. Importantly, 6 is bounded in magnitude 
by 7r/2. Although the bosonic Hamiltonian of Eq. ([7]) 
does support larger values of \S\, this parameter must 
not exceed 7r/2 in order for Eq. ([7]) to possess a fermionic 
counterpart of the form specified in Eq. ([1]). 

We proceed by converting to new charge and flavor 
fields, <I>c and corresponding to even and odd combi- 
nations of $1 and $2: 



1 



<fcix) = _[cl>i(a;) + cl>2(x)]. 
In this manner, Eqs. ([7| and ^ are rewritten as 



(10) 
(11) 



n 



i'=cj 



{V<^i,{x)fdx + edd^d 



iL|e7l[*=(«)+*/Wle-^d + H.c.} 
^|^-^[*.(o)-*,(o)]g-.^^^H.c.} 



SV2^ (dU-l]W^M 



(12) 



eV 
2^/27^ 



V$/(a;). 



(13) 



Note that the phase operator ip of Eq. ([6]) involves only 
the flavor fleld and therefore commutes with <I>c. 

At this point two successive transformations are car- 
ried out. First the canonical transformation Ti' = UTCW , 
= UYqU'' with U = exp [iV2$c(0) {d^'d- i)] is per- 
formed, to be followed by ^c{x) — > —^c{x). The latter 
operation is just a bosonic version of the particle-hole 
transformation tpi{x) "02 (^)' ^'2(2^) ^i(a^)- The 
nonequilibrium operator Yq is left unchanged by this se- 
quence of steps, and remains given by Eq. p^ . The 
Hamiltonian Ti', on the other hand, acquires one small 
but important modification as compared to Eq. (|12p : the 
parameter S is replaced with S' ^ ^ — S. As long as 
\S'\ < 7r/2, i.e., U > 0, one can revert the series of steps 
leading to Eqs. p2)l and p3)) . to recast Ti.' and Yq in 
fermionic forms. The end result of these manipulations is 
just the original pair of operators that appear in Eqs. ^ 
and ([2]), with the renormalized interaction strength 



(14) 



All other model parameters, the voltage bias included, 
are unaffected by the mapping. This proves an exact 
strong-to-weak-coupling duality of the nonequilibrium 
IRLM in the framework of bosonization.^^ 

Proceeding to observables, the duality transformation 
leaves unchanged the operators describing both the oc- 
cupation of the level, fid = d^d, and the current outgoing 
from lead j, Ij = i{e/h)y/atj[d''ipj{0) — ?/'J(0)d]. Hence 

the associated averages, Ud — {fid) and / — (Ij), obey 
the identities 

nd{T, V; Ed, h,t2,U)= nd{T, V; ea, ti,h, U'), (15) 



/(T, V; ed, ti,t2, U) = /(T, V; e^, hM, U'). 



(16) 



In particular, since U 00 maps onto [/ = 0, occupancy 
of the level and the steady-state current are identical for 
{/ ^ 00 to those of the noninteracting level. This result 
is in agreement with that of the SBA;^ including the asso- 
ciated low-energy scale being equal to F = Fi -|-r2. Here 
Fj — TTpt'j is half the tunneling rate to lead j. We empha- 
size, however, that our mapping extends the equivalence 
of [/ ^ 00 and [/ = to all moments of the current, 
implying identical full counting statistics. 

The detailed agreement with the SBA suggests that an 
analogous duality is present in the latter approach, in the 
form of TT/oC/sBA 1/{itpUsba)- Here the notation Usba 
comes to emphasize that the coupling U is defined within 
the SBA cutoff scheme. The structure of the duality is 
reflected in the expression for the low-energy scale Tq, 
which is given for irpUsBA < 1 by— 



nSBA 



{npU <l)=D[- 



t:/[tv+C{-!vpU)] 



(17) 
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with 

({x) = 2arctan(a;). 
Assuming the SBA duality then yields 



D 



7T/[7T + ai/7rpU)] 



which has the asymptotic form 



r 



l-2/(77 = p(7) 



(18) 



(19) 



(20) 



This form will be examined later in Sec. fVl 

Thus, the property of duality, Eq. (|14p . provides a sim- 
ple perspective on the nonmonotonic U dependence of 
the low-energy scale Tq, discussed in Ref. [a- Starting 
from U = and Tq — T, the low-energy scale increases 
as a function of small U r- before dropping back to its 
noninteracting value as t/ ^ oo. It must therefore de- 
crease as a function of large U, attaining a maximum 
at some intermediate coupling C/max- Assuming a sin- 
gle extremum, the maximal To must occur at the self- 
dual point [/' = U, corresponding to JJmax = 2/(7r/3) in 
bosonization and C/max — ^/{'^p) in the SBA. Indeed, 
the estimate [/max ~ 2/(7rp) is explicitly found in the 
Anderson- Yuval approachji^ Although the precise value 
of C/max may vary from one regularization scheme to an- 
other, the general estimate pC/max 1 should apply to 
all cutoff schemes, in agreement with the weak-coupling 
and numerical renormalization-group results of Borda et 
al^ Finally, we note that the derivation and structure of 
Eq. ()14|) resemble a similar duality recently found in the 
anisotropic multichannel Kondo modelJ^ 



IV. LATTICE CUTOFF 

Although the mapping of Eq. is exact within 

bosonization, it does depend on the high-energy cutoff 
scheme used. Generally speaking, all regularizations are 
expected to share the same qualitative physics. How- 
ever details, such as the parametric dependence of Tq{U), 
may vary from one realization scheme to another. It is 
quite remarkable in that respect that bosonization and 
the SBA coincide for [/ ^ oo, as both approaches em- 
ploy seemingly unrelated cutoffs. Focusing on the large- [/ 
limit, we proceed to analyze the strong-to- weak-coupling 
duality for a general lattice cutoff. Here by lattice cutoff 
we mean an arbitrary conduction-electron dispersion ek 
with a large but finite bandwidth, D. 

Consider a general lattice model for the conduction 
electrons in the leads. Using a Wilson-type construc- 
tion,^ any lattice model can be cast as a semi-infinite 
tight-binding chain with the impurity coupled to the open 
end: 



n 



j=l,2 n=0 



J = l,2 



S",oS",o 2 



(21) 



Here c, „ creates an electron in the nth Wilson shell of 
lead j . The nonequilibrium operator is given accordingly 
by 



^0 = ^ E 



„t 



(22) 



Different lattice models are distinguished by the tight- 
binding parameters e„ and which are uniquely deter- 
mined by the underlying band structure and the micro- 
scopic details of the coupling to the level. For instance, 
the case where e„ is zero for all n corresponds to particle- 
hole symmetric bands. The bandwidth D is determined 
in this case by the largest hopping matrix element ^„ 
along the chain. For convenience, we restrict attention 
hereafter to e„ = 0, taking the leads to have identical 
forms. As commented on below, both assumptions can 
be relaxed without altering the main physical result. 

Focusing onU ^ D,tj,\ed\, carry out a systematic 
expansion in 1/U. For large it is favorable energeti- 
cally to first satisfy the large capacitive coupling between 
the level c?^ and the two local lead electrons c] g with 
i — 1,2. Diagonalization of the local problem with ti 
and t2 set to zero reveals two low- lying states: c\ qc\ g|0) 
with energy — C//2, and d)\Q) with energy —U/2 + e^- 
The six remaining states are each removed in energy by 
C//2 ± e^j or more, as detailed in Table |T1 Defining the 
operator 



(23) 



it is easy to verify that (i) D projects onto the state 
d)\0); (ii) t)D^ projects onto the state c|oC2o|0); ^^'^ 
(in) i)2 = (£)t)2 ^ addition, D and D"^ anti- 

commute with all electronic operators further down the 
chain, namely, Cj^„ with n > 1. Thus, D acts as a conven- 
tional fermion annihilation operator within the truncated 
Hilbert space where only the two lowest lying local states 
are kept. 

To eliminate the locally excited states, we perform 
a canonical transformation akin to the Schrieffer- Wolff 
transformationii for the Anderson model: Ti' = e^Tie^^ . 
Here the anti-hermitian operator S is chosen such that 
the excited and low-lying sectors are decoupled within 
Ti' to linear orders in ti, t2, and ^o- The explicit form 
of the operator 5* is quite cumbersome, and will not be 
specified here. Instead, we proceed directly to the end 
result. Projecting Ti.' onto the low-energy subspace and 
settling with linear order in 1/J7, the effective low-energy 
Hamiltonian reads 



-H' 



j = l,2 n=l 



,t 
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Eigenstate 



Eigenenergy 



|0) 
dt|0) 

4,0 |o) 

4,0 |o) 
rf^4,olo) 
rf^4,olo) 

d^4,o4,olo> 



t//2 
-(7/2 + ei 





-(7/2 

f//2 + Ed 



TABLE L Eigenstates and eigenenergies of the local problem, 
defined by the Hamiltonian terms and U in Eq. (|2ip (i.e., 
the couplings ti, and ^o are all set to zero). Here |0) de- 
notes the empty state with no particles available. For large U , 
the local spectrum consists of two low-lying state: c| qCj o|0) 
with energy —(7/2, and d^|0) with energy —U/2 + €d- The six 
remaining states involve an excitation energy of [7/2 ± ed or 
more. 



where 



and 



t' = 4 



(7' = 4|^. 



(24) 



(25) 



(26) 



Here j — 3 — j marks the lead index opposite to j. To 
linear order in 1/(7, the projected nonequilibrium oper- 
ator Yq = e^Yoe~^ remains given by Eq. (|22p. apart 
from the summation over n which now runs over n > 1. 
We emphasize that Yq' does acquire corrections at order 
as does the energy in Eq. However, these 

corrections are negligibly small when U is large. 

The structure of Eq. is similar, but not identical, 
to that of Eq. ([2T|) . It differs in the form of the t'^ hopping 
terms, and in the sign of the U' interaction term. Both 
differences are conveniently accounted for by the particle- 
hole transformation 



5l 



= (-1) 



J+n+l 



(27) 



which mirrors the transformation ^dx) —^c{x) used 
in bosonization. In this manner, H' and Yq regain the 
forms of Eqs. (pij) and ([H]) modulo three differences. 
First, the summation over n now runs over n > 1. Sec- 
ond, has been replaced in Eq. (PTjl with D\ which 
couples to Cj_i instead of Cj^o (j = lj2). Third, the cou- 
plings tj and U have been replaced with t'^ and U' of 



Eqs. (|25p and ((26)) . respectively. The same differences 
carry over to the main observables of interest, namely, 
the occupancy of the level and the current outgoing from 
lead j. Explicitly, the corresponding operators transform 
according to 



/' = e' 



Ije 



= i{elh% 



'i4 



(28) 



(29) 



where we have again settled with linear order in 1/(7 
and with projection onto the low-energy sector. Thus, 
the nonequilibrium IRLM defined by Eqs. and ((^ 
has been mapped for large U onto a weakly interacting 
version of the same model. This includes the physical 
content of the operators /j and n'^ describing the current 
and occupancy of the level. 

Evidently, the nonequilibrium IRLM possess a strong- 
to-weak-coupling duality also when placed on a lattice. 
However, in contrast to bosonization, mapping of large 
to small [7 is only approximate, being controlled by the 
small parameter ^o/U <C 1, and involves a significant 
reduction of the effective tunneling matrix elements ac- 
cording to Eq. ([25]) . On the other hand, the renormalized 
interaction strength of Eq. (|26p is in good qualitative 
agreement with the predictions of bosonization, which 
follows from the fact that £,q ^ D ^ l/np for conven- 
tional lattice models. 

Our discussion thus far has focused on identical leads 
with en — along the chain. We conclude this section by 
briefly commenting on the case where either the e^'s are 
not identically zero in Eq. (|2ip (marking departure from 
particle-hole symmetry), or if the leads have different un- 
derlying band structures [i.e., e„ — *■ ej^n and ^„ —^ ^j.„ in 
Eq. (pij) ]. In the most general case, Eqs. and 
are replaced with 



and 



U' 



U 



u 



(30) 



(31) 



while td acquires the additional shift td e<i — Ci.o ~ £2,0- 
Here j = 3 — j is the lead index opposite to j. Hence, 
the capacitive coupling U' becomes lead dependent when 
the leads have different structures. Another crucial point 
pertains to the particle- hole transformation of Eq. (|27p . 
which maps each lead onto a particle-hole reflected image 
of its partner. Consequently, the structure of each lead is 
generally altered by the mapping. These modifications, 
however, do not change the main physical result, namely, 
that the nonequilibrium IRLM with large U is mapped 
onto a weakly interacting variant of the same model. 
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PHYSICAL PROPERTIES FOR LARGE U 



In case of bosonization, Eq. ((36|) takes the explicit form 



Focusing on the hmit where U is large, we next exploit 
the duality between strong and weak coupling to derive 
explicit analytic expressions for the main physical quan- 
tities of interest. These include the low-energy scale Tq, 
the differential conductance, G, and the occupancy of the 
level, Ud- 

The scale Tq is extracted from the equilibrium Hamil- 
tonian 7i, with the voltage bias set to zero. Starting from 
the dual Hamiltonian Ti.' , the latter can easily be brought 
to the canonical form considered in Ref. by converting 
to the "bonding" and "anti-bonding" combinations 



tl t2 



i^aix) = ^-^Tplix) - Jf1p2{x). 



Here 



t' = \lil + tl 



(32) 
(33) 

(34) 



is the tunneling matrix element between the level and 
the "bonding" band. No tunneling takes place between 
the level and the "anti-bonding" electrons, only capaci- 
tive coupling is left. Equations (l32 |) -([34 |) pertain to the 
dual Hamiltonian obtained within bosonization. A sim- 
ilar construction applies to the lattice version of 7i', in 
which case Eq. (p4|) is replaced with 



t' 



4^0 
U 



(35) 



[see Eq. p5|)]. Note that the conversion to "bonding" 
and "anti-bonding" modes is a standard transformation 
in equilibrium. However, it looses its usefulness when a 
finite bias is applied as the two modes couple within Yq. 

In the notations of Borda et al,^ the converted Hamil- 
tonian corresponds to N — 2, Vq — t' , and U — U' . 
Hence, one can read off the low-energy scale from the 
renormalization-group equations derived by these au- 
thors, which give 



To(p(7»l) = Z? 



l-2pU' 



(36) 



with 



r 



1 Bosonization 
(4^o/?/)^ Lattice cutoff 



(37) 



Here U' is specified by Eqs. ^1)1 and ([^ for the cutoff 
schemes used in bosonization and on a lattice, respec- 
tively. In writing Eqs. ((36)) and p7|) for a lattice cutoff, 
we have implicitly assumed that the effective bandwidth 
D and density of states p are left unchanged upon con- 
verting from Cj.o to Cj^i. 



ToipU » 1) = D ( - 



(38) 



consistent with Eq. (I20p upon the simple substitution 
l/f^SBA ^ 4/C/bosonization- This relation between the 
couplings defined within the different cutoff schemes also 
showed up in the form of the assumed duality in the SBA: 
T^pUsBA ^ ^/{t^pUsba) as compared to npU — > 4:/{TTpU) 
in bosonization. It is, however, easy to check that this 
simple relation between the coupling constants defined 
within the two schemes holds only at strong coupling, 
and does not extend to weak coupling. (The relation, in 
general, is nonanalytic, and may include singular points.) 
Note that both the SBA and bosonization expressions for 
To approach T asU oo. Indeed, fixing T and increasing 
U, Eq. ([36]) reduces asymptotically to F', which equals 
F in case of bosonization. In contrast, To{U)/To{0) 
equals {4:^0 /U)"^ for sufficiently large U in case of a lat- 
tice cutoff. The latter result agrees well with numerical 
renormalization-group calculations,— where is given 
by the hopping matrix element between the zeroth and 
first Wilson shells. 

Since U' is a marginal operator that drops to zero as 
L/ — > oo, the nonequilibrium IRLM reduces asymptot- 
ically to a noninteracting level (albeit with a strongly 
renormalized hybridization width in case of a lattice cut- 
off). The effect of U' is controlled by the dimensionless 
parameter pU' \n{D /V). When made sufficiently small, 
the differential conductance and occupancy of the level 
assume standard noninteracting forms: 



^ _ Go 



rp2 



rp2 



{ta + eVl2Y+T^ (e^ - eF/2)2 + 



(39) 



and 



nd{T, V) = - - - — arctan 

2 TT I \ Jo 

IFfl ^ fed + eV/2 

— arctan 

TT F \ To 



(40) 



with 



4FiFfl 



h (Tl + Tr)^ 



(41) 



Here we have settled for brevity with zero temperature. 
The extension to finite temperature is straightforward. 
Note that the peak conductance, Gq, is independent of 
the cutoff scheme, and is identical to its value when U = 
0. This follows from the fact that t'^/t'j^ is equal to tn/tL 
in Eq. (|25p . This result for Gq is in fact more general. It 
extends to all values of U , as can be seen from a phase- 
shift analysis. 
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VI. CONCLUSIONS 

The description of quantum impurities out of equilib- 
rium remains an outstanding theoretical challenge, with 
an urgent need for benchmark results. In this paper we 
provided such a result, by showing that the nonequi- 
librium IRLM with strong capacitive coupling is equiv- 
alent to a weakly coupled version of the same model. 
When expressed in terms of the proper low-energy scale, 
the differential conductance and occupancy of the level 
reduce for large U to standard noninteracting forms. 
This result is universal, independent of the cutoff scheme 
used. Different cutoffs are distinguished by the paramet- 
ric form of To J which is nonuniversal. In this respect 
there appears to be a qualitative difference between lat- 
tice cutoffs with a finite bandwidth, and continuum-limit 
cutoffs of the type used in bosonization and the SBA, 
where the high-energy cutoff is sent to infinity. Within 
the latter schemes. To regains its noninteracting value, 
To{U = 0) = r, as [/ ^ oo. By contrast. To decays 
to zero as To{U) ~ To{U = 0)/{pU)'^ when placed on 
a lattice. This difference must be kept in mind when- 
ever comparing different computational schemes. On the 
other hand, the peak conductance is independent of the 
cutoff scheme used, being equal to its C/ = value. 

We further wish to emphasize that the strong-to- weak- 
coupling duality reported in this paper is a distinct prop- 
erty of the IRLM with two screening channels, or leads. 
The IRLM with a single screening channel has long been 
known to map onto the anisotropic Kondo model. The 
limit U ~* oo corresponds in this mapping to — > cxd ( 
being the longitudinal Kondo coupling), while trans- 
lates to a local magnetic field. From the mapping one 
clearly sees that the single-channel IRLM with = 



flows to the same low-energy fixed point whether U is 
small or large. Namely, all values of U are physically 
equivalent at energies much smaller than To{U). How- 
ever, unlike the two-channel case, large and small U 
are inequivalent at energies comparable to and exceeding 
Tq(U). Indeed, repeating the strong-coupling expansion 
of Sec. IIVI for the case of a single screening channel one 
readily finds two low- lying local states: d^|0) and cJjO). 
As these states are occupied by a single electron each, 
charge fluctuations are strongly suppressed between the 
local complex and truncated chain at energies far below 
U. As a result, formation of a scattering center is driven 
entirely by local dynamics when U is large, in contrast 
to the case where U is small. 

A far more dramatic effect is found when the number of 
screening channels exceeds three. Here U = and U 
oo flow to different low-energy flxed points, one {U — 0) 
corresponding to a strongly hybridized level, and the 
other (U — > oo) to a free impurity.^^ The strong-coupling 
and free-impurity phases of the model are separated in 
this case by a finite-J7 Kosterlitz-Thouless transition line, 
analogous to the ferromagnetic-antiferromagnetic transi- 
tion line of the anisotropic Kondo model. A detailed 
analysis of the IRLM with multiple screening channels 
will be published elsewhere. 
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